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Abstract 

We review recent work on analytical solutions to the MSW equations for three 
neutrino flavours, for exponential and linear potentials. An application to a par- 
ticular mass matrix is also discussed. The three neutrino masses are determined, 
respectively, to be 0.001-0.004, and roughly 0.01 and 0.05 eV. 
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1 Introduction 

In this paper we shall review some re- 
cent results on analytical solutions of the 
Mikheyev-Smirnov-Wolfenstein (MSW) ef- 
fect for the propagation of three neutrino 
flavours. For two model densities which are 
relevant to the neutrino propagation in the 
sun, such results have been obtained. These 
are the exponential density and the linear 
density 0. (The case of a constant density, 
which is of some relevance for propagation 
through the Earth, has also been studied re- 
cently §.) 

For the exponential density, the solu- 
tions for the three neutrino wave functions 
H] can be expressed in terms of generalized 



hypergeometric functions, 2F2 and 3-F1. For 
the linear density, the solutions can be ex- 
pressed Fourier transform of a rather 
simple expression, which, in the case N = 
2 (two flavours) reduces to the well-known 
parabolic cylinder functions or confluent hy- 
pergeometric functions |Q. 

We also briefly discuss the application 
of the results to the Lehmann-Newton-Wu 
(LNW) mass matrix [|, 0. 



2 Exponential density 

Neutrino propagation through a medium 
where the electron neutrino (denoted 0i(r)) 
interacts differently from the others, is gov- 
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erned by the equation 
dr 
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We again introduce a new variable, z = 
ie~^. Then, the solutions to Eq. (^) can be 
expressed in terms of solutions to 
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where the mass matrix is real and symmet- 
ric. Ml = M2 = Here, D{r) = 
\/2GFNe{r), with Gp the Fermi constant 
and iVe(r) the solar electron density. 

For the sun, the density is well ap- 
proximated by an exponential, 

iVe(r) =iV,(0)e-'-/^«, ro~O.lxi?0. (2) 

It is convenient to introduce a new radial 
variable: u = r /vq + mq, and perform a rota- 
tion on the second and third components. 
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namely generalized hypergeometric func- 
tions HI: 
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Eq. (|T]) then takes the form 

. d 

du 



The 2F2 can be defined in terms of the series 
expansions 
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Let /ii, /i2 and /is be the eigenvalues of 
the 3x3 matrix 



^1 X2 X3 
X2 ^2 

Xs UJ3 



(5) 



These fij are the squares of the neutrino 
masses multiplied by ro/(2p). Together with 
uJi and UJ2 they control the evolution of the 



where {a)k is a Pochhammer symbol, {a)k = 
a{a -|- 1) ... (a -|- A; — 1). The solutions to 
Eq. (^) are thus 

V^, = Ci^«+C2^f +C3^f , (9) 

where the constants Cj are determined by 
the boundary conditions: ipi{r = 0) = 1, 
Mr = 0) = 0, Mr = 0) = 0. 

For neutrino masses and energies of phys- 
ical interest, the parameters and arguments 
of the 2F2 become too large for the series 
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expansion to be useful. A practical proce- 
dure is to adopt a stationary-phase approx- 

(3) 

imation for the 2F2 of ip\ , and express the 
others by 3F1 functions. This procedure is 
outlined below. 

Consider the ordinary differential equa- 
tion for 2F2 in the form 



For large parameters and argument, the 
Pochhammer contour V is useful [01: 
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This can be converted to an ordinary dif- 
ferential equation for 3F1 by substituting 

z = z~^: 
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This approach leads to an accurate and ef- 
ficient evaluation of the three neutrino wave 
functions in terms of stationary phase ap- 
proximations to the 2F2 and 3-F1 functions. 



3 Linear density 
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Two solutions are of the form 



The terminology "linear electron density" is 
used to mean that Ne{x) is a linear function 
/ = 0. (11) oix. 
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(12) 3.1 Two generations 



The full solutions (^ can then be con- 
structed schematically (leaving out powers) 
as follows: 



B,,F^{1) + Q,F,{2). 

(13) 

The series expansion for 3F1 has zero radius 
of convergence. However, it can be expressed 
in terms of an integral involving the familiar 
hypergeometric function 2F1, 



The case of two states and a linear potential 
has been studied extensively, starting with 
Landau and Zener in the 1930s and applied 
to neutrino mixing in the 1980s After a 
suitable scaling and shift of the variable, one 
has 
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or more explicitly: 
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-tipi{t) + a2ip2{t), 
a2Mt)- (15) 
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Elimination of ip2 gives 



-it' 



dt 



+ {al-i)iJi{t) = 0. 



(16) 



The first- derivative term can be removed 
by taking 



01 (t). 



(17) 



Then the equation for 0i(t) is 
d^Mt) 



dt^ 



+ (it' + a2-lz)0i(t) =0. (1^ 



Two hnearly independent solutions of this 
equation are the parabolic cylinder functions 



(19) 
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where p 

Parabolic cylinder functions are special 
cases of the confluent hypergeometric func- 
tion, 

(20) 

In terms of confluent hypergeometric func- 
tions \l/ and $: 
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Unfortunately, it is not clear how to gener- 
alize this approach to > 3. 

Let us therefore consider an alternative 
solution to the N = 2 case 0. We start by 
writing 
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Then it follows from Eq. (|l|) that F{() sat- 
isfies the first-order differential equation 

-C'^(C) - i:l-<nC)] + (4 - OF(C) = 
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or 
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Integrating over (, one finds 
F(C) = const. e-'«'/' 



al-e). (24) 
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With the notation [cf. Eq. (|30|)] 

hi = —00 and b^+i = +00, (26) 



the two solutions can be written as 

rbn+i 

'br 
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for n = 1,2. It can be shown that they are 
confiuent hypergeometric functions of the 
correct parameters and argument [i.e., iden- 
tical to Eq. (0)]. 



3.2 General N 

This second approach outlined above has the 
advantage that it can be generalized to an 
arbitrary number of neutrino fiavors . We 
start out by writing the equation analogous 
to (|l|) in dimensionless standard form 



z^^(t) = A(t)^(a;), 



(28) 



where 



(29) 
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and 



4 The LNW mass matrix 



A{t) 



—t a2 as ... ttN 
a2 62 ... 
aa 63 ... 



_aN . . . Bn 
There are two kinds of equations: 
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and, for A; = 2, 3, 4 ... A^: 



One finds 
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While these solutions are reasonably sim- 
ple, a numerical Fourier transform is re- 
quired. Ast ^ ±cxD, exphcit expressions can 
be written out for the different How- 
ever, these can not directly be used for im- 
posing the boundary conditions, since they 
would correspond to negative density at t ^ 



For quark mixing, it was found that a 
particular, simple texture for the d {d, s, b) 
and u {u, c, t) quark mass matrices leads to 
an acceptable CKM matrix [l^. This same 
mass matrix has been applied to the case of 
three neutrinos 0, and rather good fits to 

neu- 



the atmospheric |]TT 



and solar 
trino data have been obtained. 

The mass matrix is assumed to have the 
form 



M 



d 
c 
b 



(35) 



with 6^ = 8c^. The eigenvalues are given by 
mi, m2, and m^, with mi < m^. 

In order to outline the diagonalization, 
whereby M = RM^^i^gR^, let us introduce 
the notation 



5*1 = 171^-7712 + mi, 

= a + c 

-S2 = m-j,m2 — m^mi + m2mi, 

= 8c^ + d'^ -ac 

-S3 = mim2ms 

= ad\ 



(36) 



Then, a cubic equation for the parameter a 
can be written as 



- nSiO^ + (85^ 



52)a-53 = 0. (37) 



-00. 



A physical solution requires a real and posi- 
tive. This is equivalent to having three real 
solutions for a. One of these is negative and 
two are positive. At any point inside the al- 
lowed domain in the mi/m3-m2/m3 plane, 
there are thus two allowed solutions, denoted 
Solutions 1 and 2. 



5 



Let us consider first the atmospheric neu- 5 Summary 
trino data. The Super-Kamiokande results 
O give Am^ ^ (2 - 3) x lO'^ eV, with 



sin (26) ~ 1. The survival of muon neutri- 
nos is given by 



{31 



where U is the neutrino mixing matrix. In 
the limit of Am2it/4p <^ 1 this simplifies, 
and invoking further unitarity, one finds 



which suggests that one needs If/^sf/rsl = 
(9(1). This can be achieved within the model 
(for both Solutions 1 and 2), for mi -C m^, 
with also m2 small compared with m^. Fur- 
thermore, the data suggest that the scale 
must be such that m| ^ (2 — 3) x 10~^ eV. 

Fits to atmospheric data confirm this 
qualitative analysis. Invoking also the so- 
lar CI, Ga and Super-Kamiokande neutrino 
data 
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one finds that both Solutions 1 
and 2 give good fits for mi <^ ma, with m2 
also small as compared with m^. Forming 
a from these different atmospheric and 
solar survival probabilities, we found good 
fits [0, with m3 of the order of 0.05 eV, m2 
about 0.01 eV, and mi ~ 0.001-0.004 eV. In 
terms of the more conventional two-flavour 
analyses for the solar-neutrino sector, these 
fits roughly correspond to the large-mixing- 
angle solution. 



We have reviewed analytic work on the so- 
lutions to the MSW equations for three neu- 
trino flavours. Such results are very valuable 
for a fast scanning over the parameters of 
some given model for the mass matrix. 

Also, we have more briefly reviewed the 
LNW mass matrix, as applied to the neu- 
trino data. This is a very constrained model 
that in the quark sector describes the CKM 
matrix, and in the neutrino sector gives the 
mixing in terms of the mass eigenvalues. 

The solar neutrino data has also been 
studied within the same model, using nu- 



merical integration methods (no 2-F2's) |17 |- 
An additional fit was then found at mi ~ 
2.8 X 10~^ eV, corresponding to the small- 
mixing-angle solutions. However, this point 
is disfavoured by the atmospheric neutrino 
data, and by the electron recoil spectrum. 
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